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ABSTRACT 


A comparison of a high speed container ship using a rudder versus a podded 
propulsor is made to study replacing a rudder with a pod. A mathematical model is 
altered to simulate a ship operating with a rudder and with a pod to maneuver. The model 
incorporates the nonlinear maneuvering equations and couples the surge and sway forces, 
yaw and roll moment, and the roll angle induced during a steady turn with varying rudder 
and pod angles. The model uses the hydrodynamic derivatives and coefficients for a high 
speed container ship. The equations are numerically integrated in order to predict the roll 
angle, sway and surge velocities, and the ship’s position in the xy-plane. Both transient 
and steady state results are utilized to quantify the relative efficiency of each system. The 
results are used as a preliminary study into replacing a rudder on a ship with a podded 
propulsor. The results indicate that the ship responds faster and has a shorter turning 
radius with the pod at lower initial speeds and pod angles, while the rudder responds 
better at high speeds regardless of angle. Further research is necessary to study the effects 


of changing the pod’s position and increasing the number of pods used. 
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I. INTRODUCTION 


A. GENERAL 


Most modern ships use a rudder to provide the necessary forces and moments to 
allow the ship to turn. However, in recent years the development of podded propulsors 
has allowed ships with pods to increase their propulsive efficiency and maneuverability 
[Ref. 1]. The use of pods allows the ship to operate with a diesel-electric power plant. If 
the pods are aft on the ship, the pod can increase the forces generated while maneuvering 
at low speeds and eliminate the need for stern thrusters [Ref. 1]. Podded propulsors are 
currently used in cruise ships, where passenger comfort is the main concern, dredging 
vessels, and naval vessels such as mine hunters. The pod allows the ship to maneuver 
more rapidly with less adverse effect on the ship during turns, especially at slow speeds. 
Thus the ship can turn faster at low speeds and will have smaller roll angles. The ship 
also has a smaller turning diameter at slow speeds and pod deflections [Ref. 2]. Another 
advantage of pods is that they can be steered over 360°, allowing the ship to maneuver in 
any direction fairly easily. Pods also take up less space on a ship because the need for 


long shaft-lines, conventional rudders, and reduction gears is eliminated [Ref. 2]. 


Pods operate in a similar fashion as rudders. To turn a ship, one simply positions 
the rudder at the desired angle to accomplish the turn at a certain rate. A ship with a pod 
can apply the necessary angle by rotating one or all of the pods, if more than one pod is 
installed [Ref. 1]. Because of the weight of the pod, turning by using a pod is considered 
to cause more strain on the equipment than steering with rudders. However, rudders 
experience stall at about 35 degrees, whereas pods do not experience stall at all due to the 
fact that the force is in the same direction as the propeller thrust [Ref. 1]. Rudders are 
also used to help maintain the course of a ship when attempting to travel in a straight line. 
The ship may veer from its intended course due to ocean waves and currents. According 
to Toxopeus and Loeff, vessels with podded propulsion perform well at course keeping 


and yaw checking [Ref. 1]. 


However, pods do have some disadvantages. Pods can respond very rapidly and 


therefore may respond faster than a human operator, thus making it difficult for a 


1 


helmsman to control the ship [Ref. 2]. One solution may be to implement computer 
controls to allocate the thrusts of the propeller in such a way that the environmental loads 
do not have a detrimental effect on the ship’s performance. Pods also tend to cause large 
heel angles. According to a study performed by MARIN, most passengers tend to panic at 
heel angles of 7°, and pods have exhibited heel angles of up to 25° in some ship models 
[Ref. 2]. Also, in order for the pod to operate in the most effective manner possible, the 
ship hull has to be altered and usually more than one pod is required to provide adequate 
power and maneuverability [Ref. 1]. However, because pods have proven so versatile, 
further study should be done to see if they can be used in more applications. 


B. SCOPE OF THIS THESIS 


The purpose of this thesis is to compare the performance of a ship using a pod 
with that of a ship using a rudder with different speeds and rudder and pod deflections. 
The comparison is made to help determine the viability of replacing rudders on naval 
vessels with pods. Since naval vessels may need to maneuver rapidly at high speeds, for 
example to avoid a torpedo, it is important to determine if a pod will perform well under 


such conditions. 


The comparison of the two ships is made using a model that incorporates data for 
a high speed container ship. The model, described in Chapter II, uses the equations of 
motion for a ship to develop equations for the derivatives of several parameters. The 
original model simulates a ship with a rudder. The model is modified to simulate the use 
of a podded propulsor rather than a rudder. It is important to note that the modifications 
made on the model do not take into account any change in the hull that may be required 
for use of a pod, nor the possible need of more than one pod to properly maneuver. The 
model simply removes the rudder and assumes that the propeller, whose characteristics 


have not been altered, can now rotate in the same manner as a pod. 


The analysis performed will help determine the general operating characteristics 
of a podded propulsor over a range of speeds and pod deflections. A description of the 
transient analysis performed is found in Chapter II. The steady state analysis performed 
is described in Chapter IV. The parameters for each ship that are compared are: the total 


speed U, the sway velocity v, the roll angle , and the ship’s position in the xy-plane. 


2 


These parameters are plotted for each ship configuration and are then used to study the 


pod’s operating characteristics. 


The conclusions drawn from the analysis performed are shown in Chapter V. 


Areas of further research are also recommended. 
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I. MATHEMATICAL MODEL 


A. EQUATIONS OF MOTION 


The nonlinear equations of motion in a body-fixed coordinate system for a ship 


are shown in equations (2.1). The equations assume symmetry along the longitudinal axis 


(yc=0). The orthogonal, right-hand coordinate system is shown in Figure 1 and the 


parameters are defined in Table I. 


Surge: m| ti — Xg(q° tr’ )+z go (q+ pr)- rv-+wg |=X 


Sway: m| v+x_(pqtr)+z Z, (gr — p)+ur—wp |=Y 
Heave: m| w+ x, rp — q)-z AG +q *)+ pv-qu]=Z 


( 
I.p+(J, -1,)q bods (pq +F)- mz, (v+ur—wp)= 


Roll: % 
K-AGZ(¢) 


1,q-(1,-1,)rp-1,,(1? —p”)+ 


m| Z, (ui +qw-rv)-x, (w+ pv—qu) | =M 


Pitch: 


Yaw: Lr+(1,-1,)pqt+l,, (qr— p)+mx,(v+ur—wp)=N 


(2.1.a) 
(2.1.b) 


(2.1.c) 


(2.1.d) 


(2.1.e) 


(2.1.f) 





Figure 1. Coordinate System [Ref. 7] 


Table 1 Equations of Motion: Parameters 




































































Parameter Description 
X, y,Z Distance along the principal axes 
u, V, W Translational velocity components of ship relative to fluid 
along body axes 
p.q.r Rotational velocity components of ship relative to inertial 
reference system along body axes 
X, Y,Z Hydrodynamic force components along body axes 
K, M,N Hydrodynamic moment components along body axes 
w, 9, > w yaw angle: bow to starboard positive 
§ pitch angle: bow up positive 
droll angle: _ starboard down positive 
m Mass of ship 
XG, YG, ZG Coordinates of the center of gravity in the body axis 
system 
Ix, ly, Iz Moments of inertia about the body axis system 
Ixz, lyz, Ixy Products of inertia about the body axis system 
V Displacement volume of ship 
A Displacement weight of ship 
GZ(o) Righting moment as a function of roll angle 
6 Rudder angle in radians 
V Initial velocity of ship 
p Mass density of sea water 
L Ship length between perpendiculars (LBP) 




















B. SIMPLIFYING ASSUMPTIONS 


To simulate a ship maneuvering in the ocean without changes in depth, the 


following assumptions were made [Ref. 3]: 


1. 


The rotational velocity and acceleration about the y-axis are zero. (g =0 
and q=0) 

The translational velocity and acceleration in the z direction are zero. 
(w=0 and w=0) 


The vertical heave and pitch motions are decoupled from the horizontal 


plane motions. 


The product of inertia Ixz is small and can be neglected. 


5. The surge equation is substituted by an algebraic equation which is a 


function of u, V, and 6. 


6. The longitudinal center of gravity, (LCG), and the longitudinal center of 
buoyancy, (LCB), are at midship. 


7. The vertical center of gravity, (VCG), is on the centerline. 


8. The only important forces and moments acting on the ship induced by the 
rudder (or pod) are due to rudder (or pod) deflection. Forces and moments 


due to Sand 6 are negligible. [Ref. 4] 


Applying these assumptions to Equations (2.1), the following simplified equations 


of motion are obtained [Ref. 5]: 


Sway: m [v Pier Zap yr ur | = (2.2.a) 
Roll: 1,p—mz,(v+ur) = K —AGZ(¢) (2.2.b) 
Yaw: Li+mx,(v+ur)=N (2.2.c) 


To simplify the analysis even further, equations (2.2) are converted to a 


nondimensional form. The nondimensional relationships are shown in Table 2. 


To demonstrate the conversion to nondimensional terms, the appropriate values 


from Table 2 are substituted in Equation (2.2a) to yield the following: 
1 3 Vv? s ’ y 5 ’ Vv? . V 1 2772 
— pL |\m'|| — i'+(Lx —r' -(Zz — p' |+(Vu')| —r' ||=| ~pLV™ lY' 
[Set | (4) ae Nee Og oe 
Simplifying by factoring out (V’/L) gives 
1 2472 lon ter rey roe} 1 2472 ’ 
ae M\|V+tXor —Zop tur |= etal Y 


Equation (2.2a) can then be rewritten as Equation (2.3) as shown. 


mi] i +xgh 2g pl +ur'| Ye (2.3) 


Because Equations (2.2) and (2.3) are of the same form, the prime notation will no longer 


be included and all equations will be treated as nondimensional unless otherwise 














specified. 
Table 2. Nondimensional Parameter Relationships 
pw v=o 2 ee _ al 
y y? ~ Perv? 
L ; 
rar 2 r=r a N=aN ete 5 
V V 1/2 pLV 
L ; 1 
o2e os (F te oy 

















r 1 Pa 1 oe, — 
i= 1 [ots] is ts) AGZ(@) = s0z(0)( 5] 




















C. FORCE AND MOMENT REPRESENTATION 


Using the method described by Abkowitz and Strom-Tejsen [Ref. 4], the sway 


force, roll moment, and yaw moment can be expressed as Equations (2.4). 


Surge Force: X = f,(u,v,r, 4,6) (2.4.a) 
Sway Force: Y = f,(u,v,r,0,0,7,6,6,6) (2.4.b) 
Roll Moment: K = f,; (u,v.7,0,,7,9,6,8,6) (2.4.c) 
Yaw Moment: N = f,(u,v,1,0,0,7,9,0,6) (2.4.d) 


Using a third order Taylor expansion for f}, f2, and f3, equations (2.5) are obtained. 
f=X,Wt+X wt XV +X +XgP +X(6)+X, (2.5.a) 


fr =YvtYitY,V+V rt YF +¥,v rt Yor 
; (2.5.b) 
Yo ae Yu ay YispV + Yor os Y(6) =f Y, 


=Kyv+K Ve+K wrt+Krt+K, +Ko+K b+Ko+ 
3 v vy vyr r rrr ¢ @ @ 


(2.5.c) 
vo’ +K gro + K(6)+K, -AGZ(¢) 


2 
Kg p+ Kgs roo 
Lp=NvtNvtN,v +N, +N t NE +N, 7 + Nyb+ Not 


; : : (2.5.d) 
Nyg¥ P+ Nyy + Nigh b+ Niggt +N(6)+N, 


The coefficients in Equation (2.5) are called the hydrodynamic coefficients and represent 


the partial derivative of the force or moment with respect to the subscripted variable. For 


2 


example, Y, represents ae and Yvrr epesenie ane . The terms Yo, Ko, and No are 
Ov év\ dr? 


the sway force, roll moment, and yaw moment induced by the propeller, respectively. 
Similarly, Y(5), N(6), and K(6) represent the force and moments induced by the rudder. 
The righting moment due to the static stability of the ship is represented by the term 
AGZ(¢). 

D. RUDDER FORCE AND MOMENT REPRESENTATION 


The expressions used to determine the rudder force and moments are shown in 


Equations (2.6) and (2.7) [Ref. 5]. The parameters are defined in Table 3. 








X (6) =CpyFy sin(6) (2.6.a) 
Y(6)=—(l+a,,) Fy cos(6) (2.6.b) 
K (6) =(1+4,,)z,Fy cos(6) (2.6.c) 

N(6)=—(x_ + 4y%q ) Fy c0s(6) (2.6.d) 

19K. \(As \cse- 
Fy -( Se (4)n sin(@,) (2.7.a) 
Vp =p +V_- (2.7.b) 

SKK, 
Up =Upe, {1+ ay (2.7.c) 
UV 

-" (2.7.4) 


Up =u| (I-wp)+e{(v+x_r) +e,r+e,,7} | (2.7.e) 


Ve =VVFC 5 VA Cyl Fst V (2.7.f) 
a, =S+tan”™ [2] (2.7.2) 
Up 
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Table 3. Rudder Force and Moment Parameters 






















































































Parameter Description 
ay Rudder to hull interaction coefficient 
Fn Normal force action on the rudder 
ZR z coordinate of point on which rudder force Y35 acts 
XR x coordinate of point on which rudder force Ys acts 
XH x coordinate of point on which normal force Fy acts 
A Rudder aspect ratio 
Ar Rudder area 
VR Effective rudder inflow velocity 
OR Effective rudder inflow angle 
Ur, VR Components of rudder effective inflow velocity 
€ Constant in Equation (2.7c) 
Up Effective propeller inflow velocity 
k Constant in Equation (2.7c) 
Ky Thrust coefficient 
J Advance coefficient 
n Number of propeller revolutions per second 
D Propeller diameter 
Wp Effective propeller wake fraction 
T Constant in Equation (2.7e) 
Xp x coordinate of propeller position 
Cpvs Cpr Propeller flow rectification coefficients 
y Flow rectification coefficient 
Corrs Céprrs C8prv Rudder wake coefficients 














E. FORCES AND MOMENTS FOR CONTAINER SHIP WITH RUDDER 


The following equations are used to determine the forces and moments acting on 
the ship, including the forces and moments caused by the rudder. Equations (2.8) are used 
in the mathematical model and they combine all of the above equations with the rudder 
[Ref. 6]. The values of the hydrodynamic derivatives and coefficients are given in Table 
4. 

X=X,uw t+(1-t)T +X, wrt XV +X, 0 


Surge Force: (2.8.a) 
+X B+ Cay Py sin(5)+(m +m, )vr 
Y=Yv+Vr+Y,p+Y,O+Y,v +Y,0 +Y,,.V4r 
Sway Force: +Y,,vr? + Yuv @ + Y vo + Yo + Y lo + (2.8.b) 


(1+a,,) Fy cos(o)—(m+m, )ur 
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K=K,v+K,r+K,ptK+K,,V tK,.0 + 
Roll Moment: K,,.vr+K,,.vr° + Kv ¢ te K vb a Kg0 Se (2.8.c) 
K ul —(1+ dy )ZpFy cos(5)+m,l,ur-W (GM )¢ 
N=N,v+N,r+N,p+Nb+N, Vv +N 0 + 
Yaw Moment: Nyy Vt+N,,0° + Nygv b+ Ny gv? + (2.8.d) 
Nil P+ NygstO + (Xp t+AyXz) Fy cos(6) 
POD FORCE AND MOMENT REPRESENTATION 


The expressions used to determine the pod force and moments are shown in 


Equations (2.9). These equations were derived by comparing the thrust term in the 


equations of motion to the rudder term. The thrust terms were then modified by the same 


sign, distance (for moments), and trigonometric relationship as the rudder term, and the 


rudder term was removed from the equations. 


X, =(1-1)Tsin(6,,,) (2.9.a) 
Y, =(1-1)T cos(6,,.) (2.9.b) 
K,=—Zp (1-1)T cos(6,,,) (2.9.c) 
N, =+(I=1)T 605 (5,.) (2.9.4) 


G. FORCES AND MOMENTS FOR CONTAINER SHIP WITH POD 


The following equations are used to determine the forces and moments acting on 


the ship, including the forces and moments caused by the pod. Equations (2.10) are used 


in the mathematical model and they combine all of the above equations with the pod. The 


values of the hydrodynamic derivatives and coefficients are given in Table 4. 


X=X uw +(1-1)Tsin(5 


pod 


2 
)+X r+ Xv + 


Surge Force: (2.10.a) 
Oe + Xue” +(m+m, )vr 
ce (1-2)Tcos(6,,,)+¥.v+¥r+Y¥,p+¥,o+ 
Sway Force: Yor +Y,vrt+VY,,wr + Yuv @ + Yiugv + (2.10.b) 


Yor O+Y gro —(m+m, )ur 
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K=-z, (1-1) Tcos(5 


pod 


)+Kv+Kr+K,pt+KO+K,,v" + 
Roll Moment: Ke+K,,vr+K, vr + Kev ¢ + Kgyv + (2.10.c) 
Kr 9 + K sagt +m,l.ur-W (GM )¢ 


N=S(1-1) 1 605(3ypg)+NVAN FEN, D+ Nob 


Yaw Moment: N,vtN PP t+N wrt Nv t (2.10.d) 
Nig + NppV + Nt + N gpl 
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Table 4 


Hydrodynamic Derivatives and Coefficients [Ref. 3] 


















































































































































A. Hull only 
Xy -0.0116000 Ny -0.0038545 Ky 0.00030260 
Y; 0.0024200 N, -0.0022200 K, -0.00006300 
Yo -0.0000630 No -0.0001424 Ky -0.00002100 
een -0.1090000 Neg 0.0014920 Kwy 0.00284300 
Vee -0.0405000 Nur 0.0015600 Ker 0.00105650 
Noni 0.0017700 Nig -0.0022900 Kar -0.00004620 
ae 0.0214000 Nwvr -0.424000 Kwr -0.00055800 
Vey 0.0460500 Now -0.0190580 Ky -0.00120120 
Yy 0.0030400 Ny -0.0053766 Ky -0.00007930 
Vin 0.0093250 Nir -0.0038592 Kun -0.00024300 
x -0.0013560 N; 0.0024195 K, 0.00003569 
Xvi -0.0004226 Xvr -0.00311 Xw -0.00386 
Xu 0.00020 Xoo -0.00020 XG 0.0 
m 0.00792 mx 0.000238 my 0.000238 
I 0.0000176 l 0.0 L 0.000456 
Jx 0.0000034 Qty 0.05 Jz 0.000419 
lx 0.0313 ly 0.0313 
B. Propeller and rudder 
aH 0.237 q 1.090 (1-t) 0.825 
XH -0.480 € 0.921 XR -0.500 
Crx 0.710 k 0.631 Xp -0.526 
ZR 0.033 Car -0.156 Kr 0.527-0.455J 
Crr -0.156 Creer -0.275 CRrrv 1.960 
> 
os 0.000 eis 0.275 7 ieee S 
79.10 Fn=0.2 
Cre 0.000 Cay 1.960 Np 118.64 Fn=0.3 
158.19 Fn=0.4 
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H. 


DIMENSIONAL STATE DERIVATIVES 


Equation (2.11), which defines the derivatives of the desired variables, is used to 


find the variation of the following parameters over time: surge velocity u, sway velocity 


v, yaw velocity r, pitch rate p, roll angle , heading angle y, position in x direction, 


position in y direction, rudder deflection angle 6 (or dpoa), and n (the actual shaft 


velocity). Using Euler Integration to find the desired variable, the derivative is multiplied 


by a specified time step and then added to the previous (or initial) value of the variable. 


The variable is then equal to the new calculated value. Only u, U, v, r, , and x and y 


were used for analysis. The equation is in dimensional form. 





s &x & GH Qew BW << & 


— 





| 


-( 


(pom 


(1,45, )(1, +J,)¥ +(-m,1, )(1, +J,)K +(m,a, (1, +,)N 





—m,],)(1, + J,)¥+K(m+m, )(I,+J,)-K(m, 


(m+m, )(1,+J,)(1,+J.)-(-mJ,) (1,+5,)-(m,a,) (1,+J,) 





—(m, 


+(=myl, )(m, ane) 


(m+m, )(J, +J,)(L, +J,)—(-m, 








a,) 
ht (1,+J,) —(m, a,) (I, 

( 

| 


a,)(1,+J,)¥ +(-m,1, )(m, a ee ee ) s)= N(- =Al (7) 
(m+m,)(1,+J,)(1.+J.)-(-m, Ti ih +J,) ma ail f (1, +J,) E 
V 
it 
y 
cos(¢)r *) 


(cos(y)u —sin(y)cos(¢)v) V 
(sin (yw + cos(y) cos(¢)v) V 
-—6 


commanded 





n)x 60 


( N ommanded a 


S| 


J 


(2.11) 


The commanded rudder/pod angle is the desired deflection angle of the rudder/pod for 


the turn. The actual angle is initially zero since the ship is traveling in a straight line, and 


15 


will increase each time step as 6 is integrated. The results obtained using Euler 
Integration are then plotted to observe and compare the ship response to a rudder versus a 
podded propulsor. The parameters compared are the total speed over time, the sway 
velocity over time, the roll angle over time, and the ship’s position in the xy-plane. 
Analysis of the ship’s position in the xy-plane enables a comparison of the turning path 
of each ship configuration. The definitions for the turning path of a ship are shown in 


Figure 2. 
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Figure 2. Turning Path of a Ship [Ref. 7] 
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Hil. TRANSIENT ANALYSIS 


A. INPUT PARAMETERS 


The ship’s behavior is simulated at various initial speeds U and the rudder (or 
pod) angles as shown in Table 5. 


Table 5 Values of Varied Parameters 












































Speed U [knots] Rudder or Pod Angle 6 [deg] 
5 5 
10 10 
15 15 
20 20 
25 
30 














At each speed, the rudder and pod angles are set to the values shown above. The 
results are then plotted. The following parameters are plotted: total speed U versus time, 
roll angle » versus time, sway velocity v versus time, and position in the xy-plane. Each 
parameter is then studied at the varying speeds and rudder/pod angles. 


B. TOTAL SPEED U VERSUS TIME 


For the ship with a pod, the total speed U increases as pod angle increases (with 
the initial speed constant). The speed U does not vary as initial speed increases (with 
constant pod angle). For the ship with a rudder, the total speed U decreases as rudder 
angle increases (with constant initial speed). The speed U again does not change as initial 
speed increases (with constant rudder angle). The behavior of the speed for each ship is 
shown in Table 6 for each speed and deflection angle modeled. For both ships modeled 
the speed of the ship increases from its initial value for an initial speed of 5 knots and all 
deflections. At an initial speed of 10 and 15 knots, the ship with the pod slows down at 
deflections of 5 and 10 degrees, while the ship with the rudder speeds up (except at 15 
knots and 10 degrees where it slows down). At all other initial speeds and deflections, 


both ships will slow down. 


1, 


The model indicates that for the same initial speed U, the ship with a rudder tends 
to slow down as it maneuvers along a turn whereas the ship with a pod tends to speed up. 
This pod behavior is consistent with other studies performed on the high speed 
performance of ships with pods installed [Ref. 1]. The plots for a total speed of 5 knots 


and varying rudder/pod angles are shown in Figures 3 through 6. 
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Table 6 Ship Speed Behavior at Different Initial Speeds and 








































































































Deflection Angles 
Total Speed U Rudder or Pod Ship Speed Behavior Ship Speed 
[knots] Angle [deg] for Rudder Behavior for Pod 
5 5 il i 
5 10 1. tT 
5 15 ii t 
.) 20 tT tT 
10 e) t J 
10 10 t i} 
10 15 ii tT 
10 20 t t 
15 5 tT il 
15 10 J J 
15 15 J it 
15 20 J J 
20 5 | { 
20 10 i J 
20 15 J i) 
20 20 al A) 
2 5 J J 
2D 10 if J 
25 15 J J 
2D 20 J i} 
30 eS) i) i) 
30 10 | { 
30 15 i J 
30 20 J al 
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Figure 3. U vs. Time, U = 5 knots and 6 = 5° 
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Figure 4. U vs. Time, U = 5 knots and 6 = 10° 
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Figure 5. U vs. Time, U = 5 knots and 6 = 15° 
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Figure 6. U vs. Time, U = 5 knots and 6 = 20° 
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1600 


C. ROLL ANGLE 6 VERSUS TIME 


For the ship with a pod, the roll angle @ decreases as the pod angle 6 increases 
(with constant initial speed U). The roll angle o is approximately constant for varying 
total speed U and same pod angle 6. For the ship with a rudder, the results are the same as 
for the pod, although the magnitude of the roll angle is less than that of the ship with the 
pod for the same rudder/pod angle. The model therefore indicates that the ship with the 
pod will tend to have a larger magnitude roll angle than the ship with the rudder. This 
means that the ship with the pod will be tilted at a steeper angle during a turn than the 
ship with the rudder. These results are consistent with other studies performed on ships 
with rudders [Ref. 1]. The plots for a total speed of 5 knots and varying rudder/pod angles 
are shown in Figures 7 through 10. The plots for a rudder/pod angle of 5 degrees and 


varying total speed are shown in Figures 7 through 16. 


22 


$ [deg] 


$ [deg] 

















Roll Angle > 
U=5 knots and 6=5 deg 























2 
; : — Pad 
: : 
-10 ! 
Oo 500 1000 1500 
Time [s] 
Figure 7. o vs. Time, U = 5 knots and 6 = 5° 
Roll Angle > 
U=5 knots and 6= 10 deg 
4 
— Pod 
2D fisceesesseeeccesscsscccescesessscaesasseashesasesseesecssesssessesssesssessescssesasosyossssssscnscassssassessseqsesssesssseecse= 4 





























1 
500 1000 1600 


Time [s] 


Figure 8. o vs. Time, U = 5 knots and 6 = 10° 
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Figure 10. o vs. Time, U = 5 knots and 6 = 20° 
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o vs. Time, U = 5 knots and 6 = 5° (same as Figure 7) 
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Figure 12. o vs. Time, U = 10 knots and 6 = 5° 


25 


$ [deg] 


$ [deg] 


4 1 
a 500 1000 1600 


Rall Angle > 
U=15 knots and 6=5 deg 





— Pod 
— Rudder 














9 | L 
in} 500 1000 1500 
Time [s] 
Figure 13. o vs. Time, U = 15 knots and 6 = 5° 
Roll Angle > 
U= 20 knots and 6= Sdeg 
2 





— Pod 
—— Rudder 

















Time [s] 


Figure 14. o vs. Time, U = 20 knots and 6 = 5° 
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Figure 15. o vs. Time, U = 25 knots and 6 = 5° 
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Figure 16. o vs. Time, U = 30 knots and 6 = 5° 
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D. SWAY VELOCITY V VERSUS TIME 


For the ship with a pod, the sway velocity v is approximately constant as the 
speed U increases (with constant pod angle 5). The sway velocity v increases in 
magnitude as the pod angle 6 increases (with constant speed U). For the ship with a 
rudder, the sway velocity v is constant as the initial speed U increases (with constant 
rudder angle 5). The sway velocity v increases in magnitude as the rudder angle 6 
increases (with constant speed U). The plots for a speed of 5 knots and varying 


rudder/pod angle are shown in Figures 17 through 20. 


As seen from the plots, the sway velocity is always negative for the rudder and is 
negative for pod angles greater than 5°, regardless of the speed. However, the magnitude 


of the sway velocity is always larger for the rudder than for the pod. 
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Figure 18. v vs. time, U = 5 knots and 6 = 10° 
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Figure 19. v vs. time, U = 5 knots and 6 = 15° 
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Figure 20. v vs. time, U = 5 knots and 6 = 20° 
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E. POSITION IN XY PLANE VERSUS TIME 


The ship with a pod has a smaller turning radius than the ship with a rudder. The 
ship with the pod also has a smaller advance and transfer at the commanded pod angle for 
increasing pod angles until the total speed is 15 knots and the rudder/pod angle is 15 
degrees. After this point, the ship with the rudder has a smaller advance. Also, as the 
rudder/pod angle increases, the pod progressively has a larger turning radius than before 
for the same rudder/pod angle (and different speed U). However, the ship with the pod 
still has a smaller turning radius than the ship with a rudder. As the total speed U 
increases, the advance of the ship with the pod is greatly increased. The advance 
increases to as much as 2000 m for a total speed of 30 knots and pod angle of 20 degrees. 
The transfer also increases so that at higher speeds U the pod requires more distance to 
begin the turn. These results are consistent with previous studies [Ref. 1]. The plots for 


varying rudder/pod angle and varying total speed U are shown in Figures 21 through 32. 
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Position in XY Plane, U = 5 knots and 6 = 10° 
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Figure 24. Position in XY Plane, U = 5 knots and 6 = 20° 
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Figure 26. Position in XY Plane, U = 15 knots and 6 = 10° 
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Figure 28. Position in XY Plane, U = 15 knots and 6 = 20° 
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Figure 30. Position in XY Plane, U = 30 knots and 6 = 10° 
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IV. STEADY STATE ANALYSIS 


A. PARAMETERS OF INTEREST 


The final integrated value of each parameter found from numerical integration of 
Equation 2.11 is obtained from Matlab; this value is the steady state value of the 
parameter. For the steady state analysis, the following parameters are studied: sway 
velocity v, yaw velocity r, and roll angle 9. 


B. ANALYSIS OF PARAMETERS 


The results are plotted for different initial speeds U and different rudder/pod 
angles 5. As observed from the plots, the steady state values of each parameter do not 
vary with a change in the initial speed U (with the same rudder/pod angle). This is the 
same result obtained from the transient analysis for the parameters studied. The 
parameters vary slightly when plotted using constant initial speed U and varying 


rudder/pod angle. 
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Steady State Sway Velocity for Pod with Varying Initial Speed 


Figure 33. 
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Figure 34. 
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Figure 35. 
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Figure 38. Steady State Roll Angle for Rudder with Varying Initial Speed 
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A. 


Vv. CONCLUSIONS AND RECOMMENDATIONS 


CONCLUSIONS 


In this thesis, a mathematical model was modified to simulate the motion of a ship 


operating with a rudder and operating with a pod. The model was first used to simulate 


the ship operating with a rudder. The model was then modified to simulate operation with 


a pod. The two results were then compared to determine how well a ship with a pod 


operates at various speeds and pod deflections. Based on the comparison made, the 


following results were observed for the ship with the pod: 


1. 


The results suggest that a ship with a pod will operate very well at speeds less 
than or equal to 15 knots and pod deflections of less than or equal to 15 degrees. 
After this point, the pod tends to have a larger advance and transfer as initial 


speed and deflection angle increase. 


The results indicate that for a ship with a pod, the ship’s total speed will increase 
during a turn if the initial speed is 5 knots or less. At all other speeds and 
deflections, the ship with the pod will slow down. This behavior is similar to that 


of a ship with a rudder. 


The results show that a ship with a pod will tend to have a larger magnitude roll 
angle than a ship with a rudder; however, the roll angle does not exceed 8° 


regardless of the initial speed. 


The results indicate that the sway velocity v will be approximately constant as the 
initial speed U increases if the pod deflection is constant. The sway velocity will 


decrease as the pod deflection increases and the initial speed is held constant. 


Based on the results observed using the model as described, the rudder performs better 


than the pod at high speeds and high deflection angles. Naval vessels must have the 


capability of responding rapidly in an emergency, regardless of the initial speed of the 


vessel. The pod does not provide the necessary performance requirements to meet this 


constraint in the current configuration. 
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B. RECOMMENDATIONS 


The results produced in this thesis were intended as an initial study into the 
viability of replacing a rudder on a naval ship with a podded propulsor. In order to 
accomplish this, an existing model was modified to simulate a podded propulsor in place 
of a propeller and no rudder. The model did not take into account any change in the hull 
configuration that may be required for use of a pod, nor the possible need of more than 
one pod to properly maneuver the ship. One area of further study is to further modify the 
model used to account for this, or to create a new model to simulate a ship that already 


uses a pod. 


Further study is also necessary to determine the best possible pod and hull 
configuration for a given application. Previous studies have shown that pods do not 
operate as well as rudders at high speeds. Further study is needed to investigate the 
benefits of adding control surfaces to counteract this negative performance aspect of a 


pod. 
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APPENDIX A. MATLAB CODE FOR SHIP WITH RUDDER 


V07_CONT 


function ul = V0O7_cont(k,x,u) 


ul = [ 20*pi/180; % delta (rad) 
85]; 
V07_ DATA 


% V07_ DATA datafile for v07.m. Contains initial state 
% vector, initial input and main dimensions. 

% 

% NTH 1994 Trygve Lauvdal 

cle 

echo on 

% Default parameters for High Speed Container Ship: 
% 

% Max rudder angle : 0 deg 

% Max rudder rate _: 0 deg/s 

% Max shaft velocity : 160 rpm 


% Length of ship - 175 1 

% Initial states, inputs and sample time (to be edited by user): 
u = 15.433; % surge velocity (m/s) 
Vv = 0; % sway velocity (m/s) 
p = 0; % roll velocity (rad/s) 
r = 0; % yaw velocity (rad/s) 
phi = 0; % roll angle (rad) 
psi = 0; % yaw angle (rad) 
xpos = 0; % surge position (m) 
ypos = 0; % sway position (m) 
delta = 0; % actual rudder angle (rad) 
n = 10; % actual shaft speed (rpm) 
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dc = 20*pi/180; = % commanded rudder angle (rad) 


ne = 0; % commanded shaft speed (rpm) 

h = 0.1; % sample time (s) 

x = [u v p r phipsi xpos ypos delta n]'; % state vector 
u = [Id cnc]; % input vector 
echo 


current = pwd; 

if ~exist(['sim.m']), 

saveon = [current '/v07_data x u h']; 
else 

saveon = [current '/v07/v07_data x uh’); 
end 


eval(['save v07_data']) 
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V07_ RUDDER 


function xdot = VO7_rudder(x,ui) 

% 

% xdot = V07(x,ui); returns the time derivate of the state vector : 
% 

% X=[uvpr phi psi xpos ypos delta n ]' where 

% 


y*Uu = surge velocity (m/s) 

%vV = sway velocity (m/s) 

“Pp = roll velocity (rad/s) 

*r = yaw velocity (rad/s) 

% phi = roll angle (rad) 

% psi = yaw angle (rad) 

% Xpos = position in x-direction (m) 

% ypos = position in y-direction (m) 

% delta = actual rudder angle _— (rad) 

*n = actual shaft velocity (rpm) 

% 

% The input vector is: 

% 

wu = [ delta _cn_c ]' where 

% 

% delta _c = commanded rudder angle (rad) 
%n Cc = commanded shaft velocity (rpm) 
% 

% NTH 1994, Trygve Lauvdal 

% 

% Reference : Son og Nomoto (1982). 

% On the Coupled Motion of Steering and 

% Rolling of a High Speed Container Ship, 

% Naval Architect of Ocean Engineering, 

% 20: 73-83. From J.S.N.A. , Japan, Vol. 150, 1981. 
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% Check of input and state dimensions 
if ~(length(x) == 10), error('x-vector must have dimension 10 !");end 


if ~(length(ui) == 2),error('u-vector must have dimension 2 !');end 


% Normalization variables 


L = 175; % length of ship (m) 
sqrv = (x(1)*2 + x(2)%2); 
Vv = sqrt(sqrV); % service speed (m/s) 


% Check of service speed 


if V == 0,error('The ship must have speed greater than zero');end 


delta _max = 20; % max rudder angle (deg) 
Ddelta max = = % max rudder angle rate (deg/s) 
n_max = 160; % max shaft velocity (rpm) 


% Non-dimensional states and inputs 
delta c = ui(1); 
ui(2)/60*L/V; 


A¢ 


u =x(1IV/V;) vo =x(2)/V; 

p =x(3)*L/V;r =x(4)*L/V; 
phi = x(5); psi = x(6); 

delta = x(9);_ n =x(10)/60*L/V; 


% Parameters, hydrodynamic derivatives and main dimensions 


m=0.00792; mx =0.000238; my = 0.000238; 
Ix = 0.0000176; alphay = 0.05; Ix = 0.0313; 

ly = 0.0313; Ix = 0.0000176; Iz = 0.000456; 
Jx = 0.0000034; Jz = 0.000419; xG=0; 

B =25.40; dF = 8.00; g =9.81; 
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dA =9.00; d =8.50; nabla = 21222;%Displacement Volume 


KM =10.39; KB =4.6154; AR =33.0376;%Rudder Area 
KG =10.09; Delta = 1.8219; %Aspect Ratio 

D =6.533; %Propeller Diameter 

GM =.3/L; 

rho = 1000; t =0.175; T =0.0005; 


W~ =rho*g*nabla/(rho*L42*V%2/2); 


Xuu = = -0.0004226; Xvr =-0.00311; Xrr =0.00020; 
Xphiphi = -0.00020; Xvv =-0.00386; 


Kv = 0.0003026; Kr =-0.000063; Kp = -0.0000075; 
Kphi =-0.000021; Kvvv = 0.002843; Krrr = -0.0000462; 
Kvvr =-0.000588; Kvrr = 0.0010565; Kvvphi = -0.0012012; 


Kvphiphi = -0.0000793; Krrphi = -0.000243; Krphiphi= 0.00003569; 


Yv = -0.0116; Yr = 0.00242; Yp = 0; 

Yphi =-0.000063; Yvvv =-0.109; Yrrr = 0.00177; 
Yvvr = 0.0214; Yvrr =-0.0405; Yvvphi = 0.04605; 
Yvphiphi = 0.00304; Yrrphi= 0.009325; Yrphiphi = -0.001368; 


Nv = = -0.0038545; Nr = -0.00222; Np = = 0.000213; 

Nphi =-0.0001424;Nvvv = 0.001492; Nrrr = -0.00229; 

Nwvr =-0.0424; Nvrr = 0.00156; = Nvvphi = -0.019058; 
Nvphiphi = -0.0053766; Nrrphi = -0.0038592; Nrphiphi= 0.0024195; 


kk = 0.631; epsilon = 0.921; xR =-0.5; 
wp = 0.184; tau = 1.09; xp =-0.526; 
cpv = 0.0; cpr = 0.0; ga = 0.088; 
cRr =-0.156; cRrrr =-0.275; cRrrv = 1.96; 
cRX = 0.71; aH = _-0.237; ZR = 0.033; 
xH =-0.48; 
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% Masses and moments of inertia 
m1 1 = (m+mx); 

m22 = (mt+my); 

m32 = -my*ly; 

m42 = my*alphay; 

m33 = (Ix+Jx); 

m44 = (Iz+Jz); 


% Rudder saturation and dynamics 
if abs(delta_c) >= delta_max*pi/180, 
delta_c = sign(delta_c)*delta_max*pi/180; 


end 


delta dot = delta _c - delta; 
if abs(delta_dot) >= Ddelta_max*pi/180, 
delta_dot = sign(delta_dot)*Ddelta_max*pi/180; 


end 


% Shaft velocity saturation and dynamics 
n c=n c*V/L; 
n =n*V/L; 
if abs(n_c) >= n_max/60, 
n_c = sign(n_c)*n_max/60; 


end 


ifn > 0.3,Tm=5.65/n;else, Tm=18.83;end 
n_dot = 1/Tm*(n_c-n)*60; 


% Calculation of state derivatives 
VR =ga*v+cRr*r+ cRerr*r43 + cRrrv*12*v; 
uP =cos(v)*((1 - wp) + tau*((v + xp*r)42 + cpv*v + cpr*r)); 
J = uP*V/(n*D); 


KT =0.527 - 0.455*J; 
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uR =uP*epsilon*sqrt(1 + 8*kk*KT/(pi*J‘2)); 

alphaR = delta + atan(vR/uR); 

FN =- ((6.13*Delta)/(Delta + 2.25))*(AR/L‘2)*(uR‘%2 + vR‘2)*sin(alphaR); 
T = 2*rho*D‘%4/(V42*L“2*rho)*KT*n*abs(n); 


% Forces and moments 
X = Xuu*u’2 + (1-t)*T + Xvr*v*r + Xvv*v42 + Xrr*r*2 + Xphiphi*phir2 + ... 
(m + my)*v*r + cRX*FN*sin(delta); 


Y =Yv*v+ Yr*r+ Yp*p + Yphi*phi + Yvvv*v‘3 + Yorr*r43 + ... 


Yovr*v42*r + Yvrr*v*142 + Yvvphi*v42*phi + Yvphiphi*v*phi42 + 
Yrrphi*r’2*phi + ... 


Yrphiphi*r*phi%2 - (m + mx)*u*r + (1 + aH)*FN*cos(delta); 


K =Kv*v + Kr*r+ Kp*p + Kphi*phi + Kvvv*v‘3 + Korrr*r43 + ... 


Kyvr*v42*r + Kvyrr*v*142 + Kvvphi*v42*phi + Kvphiphi*v*phi42 + 
Krrphi*r’2*phi + ... 


Krphiphi*r*phi42 + mx*1x*u*r - W*GM*phi - (1 + H)*zR*FN*cos(delta); 


N =Nv*v+Nr*r+ Np*p + Nphi*phi + Nvvv*v43 + Nrrr*r43 + ... 


Novr*v42*r + Nvrr*v*142 + Nvvphi*v“2*phi + Nvphiphi*v*phi*2 + 
Nrrphi*r’2*phi + ... 


Nrphiphi*r*phi42 + (xR + aH*xH)*FN*cos(delta); 


% Dimensional state derivatives 


xdot =[ X*(V‘%2/L)/m11 
-(-m33*m44* Y+m32*m44*K+m42*m33*N)/(m22*m33*m44-m32%2*m44- 
m42’2*m33)*(V‘%2/L) 


(-m32*m44*¥+K *m22*m44-K*m42%2+m32*m42*N)/(m22*m33*m44- 
m32/2*m44-m42’2*m33)*(V2/L2) 


(-m42*m33*Y+m32*m42*K+N*m22*m33-N*m32’2)/(m22*m33*m44- 
m32/2*m44-m42’2*m33)*(V2/L2) 


p*(V/L) 
cos(phi)*r*(V/L) 
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(cos(psi)*u-sin(psi)*cos(phi)*v)* V 
(sin(psi)*u+cos(psi)*cos(phi)*v)*V 


delta dot 
n_dot | 

V07_SIM 
% SIM.m 
% 
% Main program for contship simulation 
% 
% NTH 1994, Trygve Lauvdal 
clear 
v07_ data % set model data 
load v0O7_ data % load model data 
N  =15000; % number of samples 
table = zeros(N,length(x)+length(u));% storage vector 
disp('Simulating'); 
for i= 1:N+1, % START OF MAIN LOOP 

table(i,:) =[x', u']; % data storage 

u =v07_cont(i*h,x,u); % input 

xdot = v0O7_rudder(x,u'); % xdot(i) = f(x(i),u(i)) 

x =x+h*xdot; % euler integration x(i+1) 
end % END OF MAIN LOOP 
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VO7RUDDER_PLOT 


function plot_data = VO7rudder_plot(table,h,N) 


figure (1) 

x = table(:,7); 

y = table(:,8); 
plot(x,y) 
legend('Pod','Rudder') 


figure(2) 

time = (0:N)*h; 

u= table(:,1)*1.94384449244; 
v = table(:,2)*1.94384449244; 
U = sqrt(u.*u + v.*v); 
plot(time,U), 
legend('Pod','Rudder') 


figure(3) 
plot(time,v), 
legend('Pod','Rudder') 


figure(4) 

phi = table(:,5)*(180/pi); 
plot(time,phi) 
legend('Pod','Rudder') 


| 
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APPENDIX B. MATLAB CODE FOR SHIP WITH POD 


V07_CONT 


function ul = V0O7_cont(k,x,u) 


ul = [ 20*pi/180; % delta_pod (rad) 
85]; 
V07_ DATA 


% V07 DATA datafile for v07.m. Contains initial state 
% vector, initial input and main dimensions. 

% 

% NTH 1994 Trygve Lauvdal 

cle 

echo on 


% Default parameters for High Speed Container Ship : 


% 
% Max pod angle : 0 deg 
% Max pod rate : 0 deg/s 


% Max shaft velocity : 160 rpm 

% Length of ship - 173 io 

% Initial states, inputs and sample time (to be edited by user): 
u = 15.433; % surge velocity (m/s) 


Vv = 0; % sway velocity (m/s) 

p = 0; % roll velocity (rad/s) 

r = 0; % yaw velocity (rad/s) 

phi = 0; % roll angle (rad) 

psi = 0; % yaw angle (rad) 

xpos = 0; % surge position (m) 

ypos = 0; % sway position (m) 

delta pod= 0; % actual rudder angle (rad) 

n = 10; % actual shaft speed (rpm) 

dc = 20*pi/180; = % commanded rudder angle (rad) 
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ne = 0; % commanded shaft speed 


h = 0.1; % sample time (s) 
x = [u v p r phipsi xpos ypos delta n]'; 

u = [Id cnc]; 

echo 


current = pwd; 

if ~exist(['sim.m']), 

saveon = [current '/v07_data x u h']; 
else 

saveon = [current '/v07/v07_data x uh’); 
end 


eval(['save v07_data']) 
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(rpm) 


% state vector 


% input vector 


V07_POD 


function xdot = V07_pod(x,u1) 

% 

% xdot = V07(x,ui); returns the time derivate of the state vector : 
% 

% X=[uvpr phi psi xpos ypos delta n ]' where 

% 


y*Uu = surge velocity (m/s) 
%vV = sway velocity (m/s) 
“Pp = roll velocity (rad/s) 
*r = yaw velocity (rad/s) 
% phi = roll angle (rad) 
% psi = yaw angle (rad) 
% Xpos = position in x-direction (m) 
% ypos = position in y-direction (m) 


% delta_pod= actual pod angle (rad) 


*n = actual shaft velocity (rpm) 

% 

% The input vector is: 

% 

wu = [ delta _cn_c ]' where 

% 

% delta _c = commanded pod angle (rad) 
%n Cc = commanded shaft velocity (rpm) 
% 

% NTH 1994, Trygve Lauvdal 

% 

% Reference : Son og Nomoto (1982). 

% On the Coupled Motion of Steering and 

% Rolling of a High Speed Container Ship, 

% Naval Architect of Ocean Engineering, 

% 20: 73-83. From J.S.N.A. , Japan, Vol. 150, 1981. 
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% Check of input and state dimensions 
if ~(length(x) == 10),error('x-vector must have dimension 10 !");end 


if ~(length(ui) == 2),error('u-vector must have dimension 2 !');end 


delta_pod = 20; % podangle (deg) 
delta pod = delta_pod*pi/180; % podangle (rad) 


% Normalization variables 


L = 175; % length of ship (m) 
sqrV = —(x(1)'2 + x(2)2); 
Vv = sqrt(sqrV); % service speed (m/s) 


% Check of service speed 


if V == 0,error('The ship must have speed greater than zero');end 


delta max = 20; % max rudder angle (deg) 
Ddelta max = > % max rudder angle rate (deg/s) 
n_max = 160; % max shaft velocity (rpm) 


% Non-dimensional states and inputs 


delta c = ui(1); 

nc = ui(2)/60*L/V; 
us = x(1)/V;; v =x(2)/V; 
p =x(3)*L/V; r =x(4)*L/V; 
phi = x(5); psi = x(6); 


delta pod = x(9); n =x(10)/60*L/V; 


% Parameters, hydrodynamic derivatives and main dimensions 


m=0.00792; mx =0.000238; my = 0.000238; 
Ix = 0.0000176; alphay = 0.05; Ix = 0.0313; 
ly = 0.0313; Ix = 0.0000176; Iz = 0.000456; 
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Jx = 0.0000034; Jz = 0.000419; xG = 0; 


B =25.40; dF = 8.00; g =9.81; 

dA =9.00; d =8.50; nabla = 21222;%Displacement Volume 
KM = 10.39; KB =4.6154; AR =33.0376;%Rudder Area 

KG =10.09; Delta = 1.8219; %Aspect Ratio 

D  =6.533; %Propeller Diameter 

GM =.3/L; 

rho = 1000; t =0.175; T =0.0005; 


W~ =rho*g*nabla/(rho*L42*V%2/2); 


Xuu = ==-0.0004226; Xvr =-0.00311; Xrr  =0.00020; 
Xphiphi = -0.00020; Xvv =-0.00386; 


Kv = 0.0003026; Kr =-0.000063; Kp = -0.0000075; 
Kphi = -0.000021; Kvvv = 0.002843; Krrr =-0.0000462; 
Kvvr =-0.000588; Kvrr = 0.0010565; Kvvphi = -0.0012012; 


Kvphiphi = -0.0000793; Krrphi = -0.000243; Krphiphi= 0.00003569; 


Yv = -0.0116; Yr = 0.00242; Yp = 0; 

Yphi =-0.000063; Yvvv =-0.109; Yrrr = 0.00177; 
Yvvr = 0.0214;  Yvrr =-0.0405; Yvvphi = 0.04605; 
Yvphiphi = 0.00304; Yrrphi= 0.009325; Yrphiphi = -0.001368; 


Nv = =-0.0038545; Nr =-0.00222; Np == 0.000213; 

Nphi =-0.0001424;Nvvv = 0.001492; Nrrr =-0.00229; 

Nvvr =-0.0424; Nvrr = 0.00156; = Nvvphi =-0.019058; 
Nvphiphi = -0.0053766; Nrrphi = -0.0038592; Nrphiphi= 0.0024195; 


kk = 0.631; epsilon = 0.921; xR =-0.5; 
wp = 0.184; tau = 1.09; xp =-0.526; 
cpv = 0.0; cpr = 0.0; ga = 0.088; 
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cRr =-0.156; cRrrr = -0.275; cRrrv = 1.96; 
cRX = 0.71; aH = _-(0..237; ZR = 0.033; 


xH =-0.48; 

% Masses and moments of inertia 
m1 1 = (m+mx); 

m22 = (mt+my); 

m32 = -my*ly; 

m42 = my*alphay; 

m33 = (Ix+Jx); 

m44 = (Iz+Jz); 


% Pod saturation and dynamics 
if abs(delta_c) >= delta_max*pi/180, 
delta_c = sign(delta_c)*delta_max*pi/180; 


end 


delta_dot = delta_c—delta_pod; 
if abs(delta_dot) >= Ddelta_max*pi/180, 
delta dot = sign(delta_dot)*Ddelta_max*pi/180; 


end 


% Shaft velocity saturation and dynamics 
n c=n_c*V/L; 
n =n*V/L; 
if abs(n_c) >= n_max/60, 
n_c = sign(n_c)*n_max/60; 


end 


ifn > 0.3,Tm=5.65/n;else, Tm=18.83;end 
n_dot = 1/Tm*(n_c-n)*60; 


% Calculation of state derivatives 


vR)  =ga*v+cRr*r + cRerr*r43 + cRrrv*r2*v; 
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uP =cos(v)*((1 - wp) + tau*((v + xp*r)*2 + cpv*v + cpr*r)); 
J = uP*V/(n*D); 

KT =0.527 - 0.455*J; 

uR =uP*epsilon*sqrt(1 + 8*kk*KT/(pi*J‘2)); 

T = 2*rho*D%4/(V42*L“2*rho)*KT*n*abs(n); 


% Forces and moments 


X = Xuu*u’2 + (1-t)*T*sin(delta_pod) + Xvr*v*r + Xvv*¥v42 + Xrr*r02 + 
Xphiphi*phi*2 + (m + my)*v*r; 


Y =(1-t)*T*cos(delta_pod) + Yv*v + Yr*r + Yp*p + Yphi*phi + Yvvv*v43 + 
Yrrr* 143 + ... 


Yovr*v42*r + Yvrr*v*142 + Yvvphi*v’2*phi + Yvphiphi*v*phi%2 + 
Yrrphi*r’2*phi + ... 


Yrphiphi*r*phi‘2 - (m + mx)*u*r; 


K =- zR*(1-t)*T*cos(delta_pod) + Kv*v + Kr*r + Kp*p + Kphi*phi + 
Kvvv*v‘3 + Korr*r43 + ... 


Kyvr*v42*r + Kvrr*v*142 + Kvvphi*v42*phi + Kvphiphi*v*phi%2 + 
Krrphi*r’2*phi + ... 


Krphiphi*r*phi2 + mx*1lx*u*r - W*GM*phi; 


N =0.5*(1-t)*T*cos(delta_pod) + Nv*v + Nr*r + Np*p + Nphi*phi + 
Nvwv*v43 + Nirr*r43 + ... 


Nvvr*v42*r + Nyrr*v*r42 + Nvvphi*v*2*phi + Nvphiphi*v*phi42 + 
Nrrphi*r’2*phi + ... 


Nrphiphi*r*phi‘%2; 


% Dimensional state derivatives 


xdot =[ X*(V42/L)/m11 
-(-m33*m44* Y+m32*m44*K+m42*m33*N)/(m22*m33*m44-m32%2*m44- 
m42’2*m33)*(V‘2/L) 


(-m32*m44* Y+K *m22*m44-K*m42%2+m32*m42*N)/(m22*m33*m44- 
m32/2*m44-m42’2*m33)*(V2/L2) 
(-m42*m33*Y+m32*m42*K+N*m22*m33-N*m32’2)/(m22*m33*m44- 
m32/2*m44-m42’2*m33)*(V2/L2) 
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p*(V/L) 
cos(phi)*r*(V/L) 
(cos(psi)*u-sin(psi)*cos(phi)*v)* V 
(sin(psi)*u+cos(psi)*cos(phi)*v)*V 


delta_dot 
n_dot ji 

V07_SIM 
% SIM.m 
% 
% Main program for contship simulation 
% 
% NTH 1994, Trygve Lauvdal 
clear 
v0O7_ data % set model data 
load v0O7_ data % load model data 
N =15000; % number of samples 
table = zeros(N,length(x)+length(u));% storage vector 
disp('Simulating'); 
for i = 1:N+1, % START OF MAIN LOOP 

table(i,:) =[x', u'] ; % data storage 

u =v07_cont(i*h,x,u); % input 

xdot = v0O7_rudder(x,u'); % xdot(i) = f(x(i),u(i)) 

x =xth*xdot; % euler integration x(i+1) 
end % END OF MAIN LOOP 


V07POD_PLOT 


function plot_data = VO7pod_plot(table,h,N) 
clf 


figure (1) 
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x = table(:,7); 

y = table(:,8); 

plot(x,y,'r') 

grid on, xlabel('X [m]'), ylabel('Y [m]') 
title('Position in the xy Plane’) 

hold on 


figure(2) 

time = (0:N)*h; 

u = table(:,1)*1.94384449244; 

v = table(:,2)*1.94384449244; 

U = sqrt(u.*u + v.*v); 

plot(time,U,'r'), 

grid on, xlabel('Time [s]'), ylabel('U [knots]') 
title(‘Total speed U') 

hold on 


figure(3) 

plot(time,yv,'r'), 

grid on, xlabel('Time [s]'), ylabel('v [knots]') 
title(‘Sway Velocity v') 

hold on 


figure(4) 

phi = table(:,5)*(180/p1); 

plot(time,phi,'r') 

grid on, xlabel('Time [s]'), ylabel(‘\phi [deg]') 
title("Roll Angle \phi') 

hold on 
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